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^ ■ ABSTRACT 

o 

\^ . We expand the off-resonant scattering theory for particle diffusion in mag- 

^ ■ netized current filaments that can be typically compared to astrophysical jets, 

including active galactic nucleus jets. In a high plasma j3 region where the di- 
rectional bulk flow is a free-energy source for establishing turbulent magnetic 
>y^ . fields via current filamentation instabilities, a novel version of quasi-linear theory 

. to describe the diffusion of test particles is proposed. The theory relies on the 

^ proviso that the injected energetic particles are not trapped in the small-scale 

O ■ structure of magnetic fields wrapping around and permeating a filament but de- 

^ . fleeted by the filaments, to open a new regime of the energy hierarchy mediated 
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by a transition compared to the particle injection. The diffusion coefficient de- 
rived from a quasi-linear type equation is applied to estimating the timescale for 
^ . the stochastic acceleration of particles by the shock wave propagating through 

the jet. The generic scalings of the achievable highest energy of an accelerated ion 



and electron, as well as of the characteristic time for conceivable energy restric- 
tions, are systematically presented. We also discuss a feasible method of verifying 
the theoretical predictions. The strong, anisotropic turbulence reflecting cosmic 
filaments might be the key to the problem of the acceleration mechanism of the 
highest energy cosmic rays exceeding 100 EeV (10^° eV), detected in recent air 
shower experiments. 

Subject headings: acceleration of particles — galaxies: jets — magnetic fields — 
methods: analytical — plasmas — shock waves 



1. INTRODUCTION 



Extremely high energy (EHE) cosmic rays beyond 100 EeV have been observed in a 
couple of decades (Takeda et al. 1998; Abbasi et al. 2004a), but their origin still remains 
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enigmatic. In regard to the generation of the EHE particles, there are two alternative 
schemes: the "top-down" scenario that hypothesizes topological defects, Z-bursts, and so on, 
and the traditional "bottom-up" (see, e.g., Olinto 2000, for a review). In the latter approach, 
we explore the candidate celestial objects operating as a cosmic-ray "Zevatron" (Blandford 
2000), namely, an accelerator boosting particle kinetic energy to ZeV (10^^ eV) ranges. By 
simply relating the celestial size to the gyroradius for the typical magnetic field strength, one 
finds that the candidates are restricted to only a few objects; these include pulsars, active 
galactic nuclei (AGNs), radio galaxy lobes, and clusters of galaxies (Hillas 1984; Ohnto 
2000). In addition, gamma-ray bursters (GRBs) are known as possible sources (Waxman 
1995). As for the transport of EHE particles from the extragalactic sources, within the GZK 
horizon (Greisen 1966; Zatsepin & Kuzmin 1966) the trajectory of the particles (particularly 
protons) ought to suffer no significant deflection due to the cosmological magnetic fleld, 
presuming its strength of the order of 0.1 nG (see, e.g., Vallec 2004, for a review). According 
to a cross-correlation study (Farrar & Biermann 1998), some super-GZK events seem to be 
well aligned with compact, radio- loud quasars. Complement arily, self-correlation study is in 
progress, showing small-scale anisotropy in the distribution of the arrival direction of EHE 
primaries (Teshima et al. 2003). More recently, the strong clustering has been confirmed, as 
is consistent with the null hypothesis of isotropically distributed arrival directions (Abbasi 
et al. 2004b). At the moment, the interpretation of these results is under active debate. 

In the bottom-up scenario, the most promising mechanism for achieving EHE is con- 
sidered to be that of diffusive shock acceleration (DSA; Lagage & Gesarsky 1983a,b; Drury 
1983), which has been substantially studied for solving the problems of particle acceleration 
in heliosphere and supernova remnant (SNR) shocks (see, e.g., Blandford &; Eichler 1987, 
for a review). In general, it calls for the shock to be accompanied by some kinds of turbu- 
lence that serve as the particle scatterers (Krymskii 1977; Bell 1978; Blandford & Ostriker 
1978). Concerning the theoretical modeling and its application to extragalactic sources such 
as AGN jets, GRBs, and so forth, it is still very important to know the actual magnetic field 
strength, configuration, and turbulent state around the shock front. At this juncture, mod- 
ern polarization measurements by using very long baseline interferometry began to reveal the 
detailed configuration of magnetic fields in extragalactic jets, for example, the quite smooth 
fields transverse to the jet axis (1803-1-784: Gabuzda 1999). Another noticeable result is that 
within the current resolution, a jet is envisaged as a bundle of at least a few filaments (e.g., 
3C84: Asada et al. 2000; 3C273: Lobanov & Zensus 2001), as were previously confirmed in 
the radio arcs near the Galactic center (GC; Yusef-Zadeh et al. 1984; Yusef-Zadeh & Morris 
1987), as well as in the well-known extragalactic jets (e.g., CygA: Perley et al. 1984, Carilli 
& Barthel 1996; M87: Owen et al. 1989). 

The morphology of filaments can be self-organized via the nonlinear development of the 
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electromagnetic current filamentation instability (CFI; Honda 2004, and references therein) 
that breaks up a uniform beam into many filaments, each carrying about net one unit 
current (Honda 2000). Fully kinetic simulations indicated that this subsequently led to the 
coalescence of the filaments, self-generating significant toroidal (transverse) components of 
magnetic fields (Honda et al. 2000a,b). As could be accommodated with this result, large- 
scale toroidal magnetic fields have recently been discovered in the GC region (Novak et 
al. 2003). Accordingly, we conjecture that a similar configuration appears in extragalactic 
objects, particularly AGN jets (Honda & Honda 2004a). It is also pointed out that the 
toroidal fields could play a remarkable role in coUimating plasma flows (Honda & Honda 
2002). Relating to this point, the AGN jets have narrow opening angles of 0oa < 10° in the 
long scales, although in close proximity to the central engine the angles tend to spread (e.g., 
(poa. ~ 60° for the M87 jet; Junor et al. 1999). Moreover, there is observational evidence that 
the internal pressures are higher than the pressures in the external medium (e.g., 4C 32.69: 
Potash & Wardlc 1980; CygA: Perley et al. 1984; M87: Owen ct al. 1989). These imply 
that the jets must be se/f coUimating and stably propagating, as could be explained by the 
kinetic theory (Honda & Honda 2002). 

In the nonlinear stage of the CFI, the magnetized filaments can often be regarded 
as strong turbulence that more strongly defiects the charged particles. When the shock 
propagation is allowed, hence, the particles are expected to be quite efficiently accelerated 
for the DSA scenario (Drury 1983; Gaisser 1990). Indeed, such a favorable environment seems 
to be well estabhshed in the AGN jets. For example, in the filamentary M87 jet, some knots 
moving toward the radio lobe exhibit the characteristics of shock discontinuity (Biretta et al. 
1983; Capetti et al. 1997), involving circumstantial evidence of in situ electron acceleration 
(Meisenheimer et al. 1996). As long as the shock accelerator operates for electrons, arbitrary 
ions will be co-accelerated, providing that the ion abundance in the jet is finite (e.g., Rawlings 
& Saunders 1991; Kotani ct al. 1996). It is, therefore, quite significant to study the feasibility 
of EHE particle production in the filamentary jets with shocks: this is just the original 
motivation for the current work. 

This paper has been prepared to show a full derivation of the diffusion coefficient for 
cosmic-ray particles scattered by the magnetized filaments. The present theory relies on a 
consensus that the kinetic energy density (ram pressure) of the bulk plasma carrying currents 
is larger than the energy density of the magnetic fields self-generated via the CFI, likely 
comparable to the thermal pressure of the bulk. That is, the flowing plasma as a reservoir 
of free energy is considered to be in a high-/? state. In a new regime in which the cosmic- 
ray particles interact off-resonantly with the magnetic turbulence having no regular field, 
the quasi-linear approximation of the kinetic transport equation is found to be consistent 
with the condition that the accelerated particles must be rather free from magnetic traps; 
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namely, the particles experience meandering motion. It follows that the diffusion anisotropy 
becomes small. Apparently, these are in contrast with the conventional quasi- linear theory 
(QLT) for small-angle resonant scattering, according to which one sets the resonance of the 
gyrating particles bound to a mean magnetic field with the weak turbulence superimposed on 
the mean field (Drury 1983; Biermann & Strittmatter 1987; Longair 1992; Honda & Honda 
2004b). It is found that there is a wide parameter range in which the resulting diffusion 
coefficient is smaller than that from a simphstic QLT in the low-/? regime. We compare 
a specified configuration of the filaments to an astrophysical jet including AGN jets and 
discuss the correct treatment of what the particle injection threshold in the present context 
could be. We then apply the derived coefficient for calculations of the DSA timescale and 
the achievable highest energy of accelerated particles in that environment. As a matter of 
convenience, we also show some generic scalings of the highest energy, taking account of the 
conceivable energy restrictions for both ions and electrons. 

In order to systematically spell out the theoretical scenario, this paper is divided into 
two major parts, consisting of the derivation of the diffusion coefficient (§ 2) and its in- 
stallation to the DSA model (§ 3). We begin, in § 2.1, with a discussion on the turbulent 
excitation mechanism due to the CFI, so as to specify a model configuration of the magne- 
tized current filaments. Then in § 2.2, we explicitly formulate the equation that describes 
particle transport in the random magnetic fluctuations. In § 2.3, the power-law spectral 
index of the magnetic fluctuations is suggested for a speciflc case. In § 2.4, we write down 
the diffusion coefficients derived from the transport equation. In § 3.1, we deal with the 
subject of particle injection, and in § 3.2, we estimate the DSA timescale, which is used to 
evaluate the maximum energy of an accelerated ion (§ 3.3) and electron (§ 3.4). Finally, § 4 
is devoted to a discussion of the feasibility and a summary. 



2. THEORY OF PARTICLE DIFFUSION IN MAGNETIC TURBULENCE 
SUSTAINED BY ANISOTROPIC CURRENT FILAMENTS 

In what follows, given the spatial conflguration of the magnetized fllaments of a bulk 
plasma jet, we derive the evolution equation for the momentum distribution function of test 
particles, which is linear to the turbulent spectral intensity, and then extract an effective 
frequency for collisionless scattering and the corresponding diffusion coefficient from the 
derived equation. 
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2.1. Model Configuration of Magnetized Current Filaments 

Respecting the macroscopic transport of energetic particles in active galaxies, there 
is direct /indirect observational evidence that they are ejected from the central core of the 
galaxies and subsequently transferred, through bipolar jets, to large-scale radio robes in 
which the kinematic energy considerably dissipates (e.g., Biretta et al. 1995; Tashiro & Isobe 
2004). In this picture, it is expected that the directional plasma flows will favorably induce 
huge currents in various aspects of the above transport process (e.g., Appl & Camenzind 
1992; Conway et al. 1993; an analogous situation also seems to appear in GRB jets: e.g., 
Lyutikov & Blandford 2003). Because of perfect conductivity in fully ionized plasmas, hot 
currents driven in, e.g., the central engine prefer being quickly compensated by plasma 
return currents. This creates a pattern of the counterstreaming currents that is unstable 
for the electromagnetic CFI. As is well known, the pattern is also unstable to electrostatic 
disturbances with the propagation vectors parallel to the streaming direction, but in the 
present work, we ehminate the longitudinal modes, so as to isolate the transverse CFI. 
Providing a simple case in which the two uniform currents are carried by electrons, the 
mechanism of magnetic field amphfication due to the CFI is explained as follows. When 
the compensation of the counterpropagating electron currents is disturbed in the transverse 
direction, magnetic repulsion between the two currents reinforces the initial disturbance. As 
a consequence, a larger and larger magnetic field is produced as time increases. For the 
Weibel instability as an example, the unstable mode is the purely growing mode without 
oscillations (Honda 2004), so the temporal variation of magnetic fields is expected to be 
markedly slow in the saturation regime (more on these is given in §§ 2.2 and 2.3). Note that 
a similar pattern of quasi-static magnetic fields can be also established during the collision 
of electron-positron plasmas (Kazimura et al. 1998; Silva et al. 2003) and in a shock front 
propagating through an ambient plasma with/without initial magnetic fields (Nishikawa et 
al. 2003). These dynamics might be involved in the organization of the knotlike features in 
the Fanaroff-Riley (FR) type I radio jets, which appear to be a shock caused by high- velocity 
material overtaking slower material (e.g., Biretta et al. 1983). Similarly, the cumulative 
impingement could also take place around the hot spots of the FR type II sources, which 
arguably reflect the termination shocks. In fact, the fllamentary structure has been observed 
in the hot spot region of a FR II source (Perley et al. 1984; Carilli & Barthel 1996). 

Furthermore, estimating the energy budget in many radio lobes implies that the ram 
pressure of such current-carrying jets is much larger than the energy density of the magnetic 
flelds (Tashiro & Isobe 2004); namely, the jet bulk can be regarded as a huge reservoir of free 
energy. In this regime, the ballistic motion is unhkely to be affected by the self-generated 
magnetic fields, as is actually seen in the hnear feature of jets. As a matter of fact, the 
GC region is known to arrange numerous linear filaments, including nonthermal filaments 
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(Yusef-Zadeh et al. 2004). 

Taking these into consideration, we give a simple model of the corresponding current- 
magnetic field system and attempt to unambiguously distinguish the present system from 
the one that appears in the low-/3 plasmas hitherto well studied. In Figure 1, for a given 
coordinate, we depict the configuration of the linear current filaments and turbulent mag- 
netic fields of the bulk plasma. Recalling that magnetic field perturbations develop in the 
direction transverse to the initial currents (e.g., Honda 2004), one supposes the magnetic 
fields developed in the nonlinear phase to be B = (B^, By, 0) (Montgomery & Liu 1979; 
Medvedev & Loeb 1999), such that the vectors of zeroth-order current density point in the 
directions parallel and antiparallel to the ^-direction, i.e., J ~ Jz, where the scalar J (^ 0) 
is nonuniformly distributed on the transverse x-y plane, while uniformly distributed in the 
z-direction. Note that for the fiuctuating magnetic field vectors, we have used the simple 
character (B) without any additional symbol such as "5," since the establishment of no 
significant regular component is expected, and simultaneously, J (~ V x B) well embodies 
the quasi-static current filaments in the zeroth order. For convenience, hereafter, the nota- 
tions "parallel" (||) and "transverse" (_L) are referred to as the directions with respect to 
the linear current filaments aligned in the they are well defined reasonably (n.b. 

in the review of § 2.4, ||{, and _Lb with the subscript "6" refer to a mean magnetic field line). 
It is mentioned that the greatly fluctuating transverse flelds could be reproduced by some 
numerical simulations (e.g., Lee & Lampe 1973; Nishikawa et al. 2003). In an actual fila- 
mentary jet, a significant reduction of polarization has been found in the center, which could 
be ascribed to the cancellation of the small-scale structure of magnetic fields (Capetti et al. 
1997), compatible with the present model configuration. In addition, there is strong evidence 
that random fields accompany GRB jets (e.g., Greiner et al. 2003). The arguments expanded 
below highlight the transport properties of test particles in such a bulk environment, that 
is, in a forest of magnetized current filaments. 



2.2. The Quasi-linear Type Equation for Cosmic-Ray Transport 

We are particularly concerned with the stochastic diffusion of the energetic test particles 
injected into the magnetized current filaments (for the injection problem, see the discussion in 
§3.1). As a rule, the Vlasov equation is appropriate for describing the collisionless transport 
of the relativistic particles in the turbulent magnetic field, B(r,t), where r = {x,y) and the 
slow temporal variation has been taken into consideration. Transverse electrostatic fields 
are ignored, since they preferentially attenuate over long timescales, e.g., in the propagation 
time of jets (see § 2.3). The temporal evolution of the momentum distribution function for 
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the test particles, /p, can then be described as 



Dt 



-f 



dr ) '^^ ^ c 



(vxB).^ 



/p 







for arbitrary particles. Here q is the particle charge/ c is the speed of light, and the other 
notations are standard. We decompose the total distribution function into the averaged and 
fluctuating part, /p = (/p) +Sfp, and consider the specific case in which from a macroscopic 
point of view, the vector B = {B^, By, 0) is randomly distributed on the transverse x-y 
plane (Montgomery & Liu 1979; Medvedev & Loeb 1999). Taking the ensemble average of 
equation (1), (D/p/Di) — 0, then yields 



(vxB).| 



(2) 



where we have used (B) ~ 0. Taking account of no mean field implies that we do not invoke 
the gyration and guiding center motion of the particles. Subtracting equation (2) from 
equation (1) and picking up the term linear in fluctuations, viz., employing the conventional 
quasi-linear approximation, we obtain 



As usual, equation (3) is vahd for (/p) » |5/p| (Landau & Lifshitz 1981). As shown in 
§ 3.1, this condition turns out to be consistent with the aforementioned implication that 
the injected test particles must be free from the small-scale magnetic traps embedded in the 
bulk. Relating to this, note that to remove the ambiguity of terminologies, the injected, 
energetic test particles obeying /p are just compared to the cosmic rays that are shown 
below to be diffusively accelerated owing to the present scenario. Within the framework of 
the test particle approximation, the back-reaction of the slow spatiotemporal change of (/p) 
to the modulation of B (sustained by the bulk) is ignored, in contrast to the case for SNR 
environments, where such effects often become nonnegligible (e.g.. Bell 2004). 

In general, the vector potential conforms to B = V x A and V ■ A = 0. For the stan- 
dard, plane wave approximation, we carry out the Fourier transformation of the fluctuating 
components for time and the transverse plane: 



fv X B] 



d_ 
dp 



(/p) 



(3) 



^For example, q = — \e\ for electrons, q = \e\ for positrons, and q = Z\e\ for ions or nuclei, where e and Z 
are the elementary charge and the charge number, respectively. 
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5/p(r,i) = J 5/p,;ce^'^^ "-)-'^*ld3/C, (4) 
A(r,i) = y"A;ce'^('*-"^"''*Jd% (5) 
B(r,i) = y"B;ce'K''"^-'"*ld2/C, (6) 

and B/c = x Ajc, where i = /C = {k, cj}, and d^/C = d^kdcj. The given magnetic 

field configuration follows A = Az (A^; = A/cz) and k ± z. As illustrated in Figure 1, the 
scalar quantity A 0) is also random on the transverse plane, with no mean value. Making 
use of equations (4)-(6), equation (3) can be transformed into 

SU,K = f [u; - (k • v)]-^ [v X (k X A^)] • (7) 
On the other hand, the right-hand side (RHS) of equation (2) can be written as 

RHS = -^^ ^1 d^/C'e^K"^'-)-'*! |[v X (k' x A,c')] ■ ^} . (8) 

Substituting equation (4) (involving eq. [7]) into equation (8), equation (2) can be expressed 
as 



d(/p/ ,y / / ^3^^3j^^/gi{[(k+k')-r]-(a;W)t} 



dt 



4{f 



/, / . 9 f [v X (k X Aic)] d } , ^ ,\ 
lvx(k'xA,01.g^|L-A_J._|(;„)^, (9) 

where the definition of the total derivative, d/dt = d/dt + v ■ (d/dr), has been introduced. 
As for the integrand of equation (9), it may be instructive to write down the vector identity 
of 



V X (kW X Aicio) = (v • Aicio) kW - (k^ • v) A^,). (10) 

From the general expression of equation (9), we derive an effective collision frequency that 
stems from fiuctuating field-particle interaction, as shown below. 
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For convenience, we decompose the collision integral (RHS of eq. [9]) including the scalar 
products, • (d/dp), into the four parts: 



dt 

where i,j = 1, 2 and 



(11) 



(|,3^a3/c'... A... A (/„)). (12) 

In the following notations, the subscripts "1" and "2" indicate the parallel (||) and per- 
pendicular (_L) direction to the current filaments, respectively. Below, as an example, we 
investigate the contribution from the integral lu (see Appendix for calculation of the other 
components). For the purely parallel diffusion involving the partial derivative of only d/dp\\, 
the first term of the RHS of equation (10) does not make a contribution to equation (9). 

In the ordinary case in which the random fluctuations are stationary and homogeneous, 
the correlation function has its sharp peak at a; = —lo' and k = — k' (Tsytovich & ter Haar 
1995), that is. 



{AicAic') = \A\lj{k + k')6{u; + u;'), 



(13) 



where the Dirac 5-function has been used. Here note the relation of l^lk^) — l^l-k -oj because 
we have 74_k,-c<; = ^k<^, where the superscript asterisk indicates the complex conjugate; this 
is valid as far as A(r, t) is real, i.e., B(r, t) is observable. By using equation (13), the integral 
component lu can be expressed as 



k- V 



d 



2 p f) 

7n = i% / d^kda;|^|^ (k-v) — 

J ' dp\\ [a; - (k ■ v) dp\\ 



(/p) 



(14) 



where the relation of A^; ■ (d/dp) = Ajc {d/dp\\) has been used. In order to handle the reso- 
nant denominator of equation (14), we introduce the causality principle o{\im.^_^Q [a; — (k ■ v) -f ie] 
P [cu — (k ■ v)]~^ — in 6 [cj — (k ■ v)], where V indicates the principal value (Landau & Lif- 
shitz 1981). One can readily confirm that the real part of the resonant denominator does 
not contribute to the integration. Thus, we have 



7rg 

^2 



d^kd^l^l^^^ 



(k.v)^{.[.-(k.v)](k.v)^(/,)}. 



(15) 



-10- 



Equation (15) shows the generahzed form of the quasi-hnear equation, allowing \A\l^^^ to be 
arbitrary functions of k and u;.^ In the present circumstances, a typical unstable mode of 
the CFI is the purely growing Weibcl mode with = in collisionless regimes, although in 
a dissipative regime the dephasing modes with a finite but small value of c<j = ±Ac<Jk are 
possibly excited (Honda 2004). In the latter case, the spectral lines will be broadened in the 
nonlinear phase. Nevertheless, it is assumed that the spectrum still retains the peaks around 
cu ~ ±Aa;, accompanied by their small broadening of the same order, where | Aci;| <^ 7k ~ i-^p, 
and 7k and cUp/ (27r) are the growth rate and the plasma frequency, respectively. In the special 
case reflecting the purely growing mode, the spectrum retains a narrow peak at a; = with 
|Aa;| ~ (Montgomery & Liu 1979). Apparently, the assumed quasi-static properties are 
in accordance with the results of the fully kinetic simulations (Kazimura et al. 1998; Honda 
et al. 2000a), except for a peculiar temporal property of the rapid coalescence of filaments. 
Accordingly, here we employ an approximate expression of 

l^lk,., ~ \A\lS {cu - Acu) + \AWS {cu + Acu) , (16) 

where \A\'^ = |A|^k- Note that when taking the limit of |Aa;| — >• 0, equation (16) degenerates 
into \A\l^r^2\A\i5i^). 

Substituting equation (16) into equation (15) yields 




Furthermore, we postulate that the turbulence is isotropic on the transverse plane, though 
still, of course, allowing anisotropy of the vectors A parallel to the z-axis. Equation (17) 
can be then cast to 



2V Q i-^g 

hi ~ 



o v±Tr- / T^cos^^ / dk27rk6 (Au - kv^ cose) k^\A\lv^— (U , (18) 
dp\i J 27: J " dp\\ 

where k — \k.\, v± — |v_l|, and k • v = k • (vj_ + v\\zj — kv± cos9. 



^For the case in which the unstable mode is a wave mode with Wk ^ 0, the frequency dependence of the 
correlation function can be summarized in the form of |JF|^ ^ = l-^lk*^ ^ ^k) + l-^l-k*^ + ^k), which is 
valid for weak turbulence concomitant with a scalar or vector potential JF. However, this is not the case 
considered here. The free-energy source that drives instability is now current flows; thereby, unstable modes 
without oscillation (or with quite slow oscillation) can be excited. 
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As concerns the integration for d, we see that the contribution from the marginal region 
of the smaller | cos^|, reflecting narrower pitch angle, is negligible. In astrophysical jets, the 
pitch angle distribution for energetic particles still remains unresolved, although the distribu- 
tion itself is presumably unimportant. Hence, at the moment it may be adequate to simply 
take an angular average, considering, for heuristic purposes, the contribution from the range 
of I cos^l ~ 0(1) ^ e for a small value of e. If one can choose e > \^uj\/ {k^i^v_i_), the above 
relation, e <C | cos 6*1, reflects the off-resonant interaction, i.e., |Aci;| <^ |k ■ v|. The minimum 
wavenumber, kraim is typically of the order of the reciprocal of the finite system size, which 
is, in the present circumstances, larger than the skin depth c/uj^. These ensure the aforemen- 
tioned relation of | Aa;| <^ uj^ (or | Aa;| ~ 0). In addition, the off-resonance condition provides 
an approximate expression of (5 (Ao; — kv±cos9) ~ [kv^] cos^^|)^^. Using this expression, the 
integral for the angular average can be approximated by ~ cos^ ^^/(27r| cos^^|)d^^ = 2/7r. 
This is feasible, on account of the negligible contribution from the angle of | cos^| < e. Then 
equation (18) reduces to 



where we have defined the modal energy density (spectral intensity) of the quasi-static turbu- 
lence by Ik = Stt/c (A;^|>l|^/47r), such that the magnetic energy density in the plasma medium 
can be evaluated hy ^ I St: ^ if"^ hdk. 



The energy density of the quasi-static magnetic fields, Um, likely becomes comparable 
to the thermal pressure of the filaments (Honda et al. 2000a,b; Honda & Honda 2002). 
When exhibiting such a higher Um level, the bulk plasma state may be regarded as the 
strong turbulence; but recall that in the nonlinear CFI, the frequency spectrum with a sharp 
peak at cl' = is scarcely smoothed out, since significant mode-mode energy exchanges are 
unexpected. This feature is in contrast to the ordinary magnetohydrodynamic (MHD) and 
electrostatic turbulence, in which a larger energy density of fluctuating fields would involve 
modal energy transfer. One of the most remarkable points is that as long as the vahdity 
condition of the quasi-linear approximation, (/p) ^ \Sfp\, is satisfied (for details, see § 3.1), 
the present off-resonant scattering theory covers even the strong turbulence regime. That is, 
the theory, which might be classified into an extended version of the QLT, does not explicitly 
restrict the magnetic turbulence to be weak (for instruction, Tsytovich & ter Haar [1995] 
have considered a generalization of the quasi-linear equation in regard to its application to 




(19) 



2.3. Spectral Intensity of the Transverse Magnetic Fields 
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strong electrostatic turbulence) . Apparently, this is also in contrast to the conventional QLT 
for small-angle resonant scattering, which invokes a mean magnetic field (well defined only 
for the case in which the turbulence is weak) in ordinary low- (3 plasmas. 

In any case, in equation (19) we specify the spectral intensity of the random magnetic 
fields, which are established via the aforementioned mechanism of the electromagnetic CFI. 
The closely related analysis in the nonlinear regime was first performed by Montgomery & 
Liu (1979), for a simple case in which two countcrstreaming electron currents compensate 
for a uniform, immobile ion background. In the static limit of ^ 0, they have derived 
the modal energy densities of fluctuating electrostatic and magnetic fields, by using statis- 
tical mechanical techniques. They predicted the accumulation of magnetic energy at long 
wavelengths, consistent with the corresponding numerical simulation (Lee & Lampe 1973). 
It was also shown that at long wavelengths, the energy density of a transverse electrostatic 
field was comparable to the thermal energy density. However, when allowing ion motions, 
such an electrostatic field is found to attenuate significantly, resulting in equipartition of the 
energy into magnetic and thermal components (Honda et al. 2000a,b). That is why we have 
neglected the electrostatic field in equation (1). 

When the spectral intensity of the magnetic fiuctuations can be represented by a power- 
law distribution of the form 

Ik oc A;-", (20) 

we refer to a as spectral index. Montgomery & Liu (1979) found that for the transverse 
magnetic fields accompanying anisotropic current filaments, the spectral index could be 
approximated by a ~ 2 in a wide range of k, that is, 

Ik oc A;"^ (21) 

Note that the spectral index is somewhat larger than q;mhd — 1 — 5/3 for the classical MHD 
context (Kohnogorov 1941; Bohm 1949; Kraichnan 1965). The larger index is rather consis- 
tent with the observed trends of softening of filamentary turbulent spectra in extragalactic 
jets (e.g., a ~ 2.6 in Cyg A; CariUi & Barthel 1996, and references therein). Although 
the turbulent dissipation actually involves the truncation of Ik in the short-wavelength re- 
gions, we simply take k^ax oo, excluding the complication. Using equation (20) and the 
expression of the magnetic energy density Wm, we find the relation of 
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rp.^^^m (22) 

for q; > 1. The spectral details for individual jets (such as the bend-over scales of Ik, 
correlation length, and so on; e.g., for the hehosphere, see Zank et al. 1998, 2004) will render 
the integration of equation (22) more precise, but the related observational information has 
been poorly updated thus far. For the present purpose, we simply use equation (22), setting 
^rnin = 7r/i?, whcre R stands for the radius of the jet, which is actually associated with the 
radius of a bundle of filaments of various smaller radial sizes (e.g., Owen et al. 1989). This 
ensures that the coherence length of the fluctuating force, ~ k~^, is small compared with a 
characteristic system size, i.e., the transverse size, as is analogous to the restriction for use 
of the conventional QLT. 



2.4. The Diffusion Coefficients 

In order to evaluate the diffusion coefficients of test particles, one needs to specify 
the momentum distribution function, (/p), in equation (19). As is theoretically known, 
the Fermi acceleration mechanisms lead to the differential spectrum of dn/dE oc 
[or n{> E) (X E~^^^; Gaisser 1990], where dn{E) defines the density of particles with 
kinetic energy between E and E + dE. For the first-order Fermi mechanism involving 
nonrelativistic shock with its compression ratio of r < 4, the power-law index reads /3 — 
(r + 2) / (r — 1) > 2, accommodated by the observational results. With reference to these, 
we have the momentum distribution function of (/p) oc |p|^*^^+^) for the ultrarelativis- 
tic particles having E = |p|c, such that in the isotropic case, the differential quantity 
(/p) |Ppd|p|/(27r^) corresponds to dn{E) defined above (e.g., Blandford & Ostriker 1978). 
Then, in equation (19) the partial derivative of the distribution function can be estimated 
as 92(/p)/9pf ~ (/3 + 2)p + 3)(p||/|p|)2 - {pJ\p\f]{cyE'){f^), where we have used 
|Pp — p| Making use of this expression and equation (22), equation (19) can be ar- 

ranged in the form of In ~ un (/p). Here i/n refiects an effective collision frequency related 
to the purely parallel diffusion in momentum space, to give 

2{a-l){(3 + 2) [{(3 + 3) - ^j] j;, cq'B'R 
^ 

where we have used the definitions of B"^ = (|Bp), and '0i = P||/|p| ^ and ■02 = P±/|p| > 0, 
whereby ^-i/jf = 1. 
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Similarly, one can calculate the other components of the integral as outlined in the 
Appendix and arrange them in the form of lij ~ i/ij (/p). As a result, we obtained 

2 (g - 1) (/? + 2) {(3 + 4) ^1^2 cq'B'R 
^ ^^^^ 

and 

1^12 = -1^11, 

1^21 = -1^22- (25) 



As would be expected, we confirm a trivial relation of d (/p) /dt = ■ lij ~ Yli j ^ij (/p) — 0) 
stemming from the orthogonality in the RHS of equation (2). 

Now we estimate the spatial diffusion coefficients in an ad hoc manner: ~ c^ipiipj/ (2z/jj). 
It is then found that the off-diagonal components, K12 and K21, include the factor of sgn(-?/'i ^ 
0) = ±1, implying that these components vanish for an average. For ipf = ^ and ipl = \ re- 
flecting the momentum isotropy, the diffusion coefficients can be summarized in the following 
tensor form: 



f K\\ \ -\/67rQ; cE'^ 

V ) ^ 8 (a - 1) q^B^R 

where n\\ = Kn and k± = K22- The perpendicular component can be expressed 

where k = 2(/9 + -|- 4). Note the allowable range of 1 < k < 2 for > 2; particularly, 

K ~ 1 for the expected range of /? ~ 2 — 3. 

It may be instructive to compare the diffusion coefficient of equation (26) with that 
derived from the previously suggested theories including the QLT. In weakly turbulent low- 
/? plasmas, the mean magnetic field with its strength B, which can bind charged particles and 
assign the gyroradius of Vg = E/{\q\B), provides a well-defined direction along the field line; 
therefore, in the following discussion, we refer, for convenience, to and ±5 as the parallel 
and perpendicular directions to the mean magnetic field, respectively. For a simplistic QLT, 
one sets an ideal environment in which the turbulent Alfven waves propagating along the 
mean field line resonantly scatter the bound particles, when k'^^ ~ Vg, where /cy^ is the parallel 
wavenumber (Drury 1983; Longair 1992). Assuming that the particles interact with the 
waves in the inertial range of the turbulent spectrum with its index af,, the parallel diffusion 
coefficient could be estimated as (Biermann & Strittmatter 1987; Miicke & Protheroe 2001) 



(/3+l)(/9+2) 





(26) 



(/3+2)(^+4) 
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1 CT 

for ctft 7^ 1 and < k^^^^^, where k'^^^^^ reflects the correlation length of the turbulence 
and rjb (< 1) defines the energy density ratio of the turbulent/mean field. In the special 
case of ctb = 1, referred to as the Bohm diffusion limit (Bohm 1949), one gets the ordering 
K\\^ ~ K.B/rjb, where kb = cVg/S denotes the Bohm diffusion coefficient for ultrarelativistic 
particles. Considering the energy accumulation range of smaller k\\^ for the Kolmogorov 
turbulence with aij = 5/3, Zank et al. (1998) derived a modified coefficient that recovered 
the scaling of equation (27) in the region of rg -C k'^^^^^. As for the more complicated 
perpendicular diffusion, a phenomenological hard-sphere scattering form of the coefficient is 
>^±b ~ Vb'^Wb Bohm diffusion limit; and Jokipii (1987) suggested a somewhat extended 

version, k±^^ — + ('^||(,/'"g)^] (referred to as kj below), where Ay^ is the parallel mean 

free path (mfp). A significantly improved theory of perpendicular diffusion has recently 
been proposed by Matthaeus et al. (2003), including nonlincarity incorporated with the two- 
dimensional wavevector whereupon for cif, = 5/3, Zank et al. (2004) have derived an 
approximate expression of the corresponding diffusion coefficient, although it still exhibits a 
somewhat complicated form (referred to as Kz)- 

On the other hand, within the present framework the gyroradius of the injected energetic 
particles cannot be well defined, because of | (B) | ~ (§§ 2.1 and 2.2). Nonetheless, in order 
to make a fair comparison with the order of the components of equation (26), the variable 
B is formally equated with i? = (|Bp)^''^. In addition, the correlation length is chosen as 
k\\^,mm ^ R~^) corresponding to the setting in § 2.3. Then the ratio of kh for k = 1 to 
equation (27) is found to take a value in the range of 



K , / 1 E 1 mGlOO pc\'''' 



Z 100 EeV B R J 



for the expected values of «, /3 ^ 2 — 3. Here k = Ku and q = Z\e\ have been introduced. 
Similarly, we get the scaling of h/hb ~ 10~^{E/ZeBR), and k/kj ~ [E / ZeBRy/^ for 
= 5/3 and r/ft ~ 10^^ followed by Ay^ ^ rg. Furthermore, considering the parameters 
given in Zank et al. (2004), which can be accommodated with the above rib ~ 10~^, we 
also have k/kz ~ {E/ZeBRy^/^ in the leading order of kz, proportional to r^J^ . These 
scahngs are valid for arbitrary species of charged particles; for instance, setting Z — 1 reflects 
electrons, positrons, or protons (sec footnote 1). Particularly, for k < k\\^ in equation (28), 
the efficiency of the present DSA is expected to be higher than that of the conventional 
one based on the simplistic QLT invoking parallel diffusion (Biermann & Strittmatter 1987). 
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This can likely be accomplished for high-Z particles, as well as electrons with lower maximum 
energies. Here note that k cannot take an unlimitedly smaller value with decreasing E, since 
the effects of cold particle trapping in the local magnetic fields make the approximation of 
no guide field (eq. [2]) worse; and the lower limit of E is relevant to the injection condition 
called for the present DSA. More on these is given in § 3. 

To apply the DSA model, one needs the effective diffusion coefficient for the direction 
normal to the shock front, referred to as the shock-normal direction. For convenience, here 
we write down the coefficient for the general case in which the current filaments are inclined 
by an angle of with respect to the shock-normal direction. In the tensor transformation of 



k' = T , (29) 

\ '^21 '^22 / 

^^/cos0 -sin0\ 
\ sm0 COS0 / 

we identify the shock-normal component with k'^. It can be expressed as 



or simply for 1, where the subscripts = I, II indicate the upstream and 

downstream regions, respectively. The expression of equation (31) appears to be the same as 
equation (4) in Jokipii (1987). However, note again that now || and ± refer to the direction 
of the linear current filaments, compared to an astrophysical jet (§ 2.1 and Fig. 1). 



3. PARTICLE ACCELERATION BY SHOCK IN MAGNETIZED 

CURRENT FILAMENTS 

We consider the particle injection mechanism that makes the present DSA scenario 
feasible, retaining the validity of the quasi-linear approximation. Then, using the diffusion 
coefficient (cq. [26]), we estimate the DSA timescale for arbitrary species of charged particles 
and calculate, by taking the competitive energy loss processes into account, at the achievable 
highest energies of the particles in astrophysical filaments. 



-17- 



3.1. The Conception of Energy Hierarchy, Transition, and Injection of 

Cosmic-Ray Particles 

In the usual DSA context, equation (31) that calls equation (26) determines the cycle 
time for one back-and-forth of cosmic-ray particles across the shock front, which is used 
below for evaluation of the mean acceleration time (§ 3.2; Gaisser 1990). Here we note that 
equation (26) is valid for a high-energy regime in which the test particles with E are unbound 
to the local magnetic fields, so as to experience the nongyrating motion. As shown below, this 
limitation can be deduced from the validity condition of the quasi-linear approximation that 
has been employed in § 2.2. Using equations (4) and (7), the validity condition (/p) ^ \Sfp\ 
can be rewritten as 



(/p) » 



^ / d^ke^"^-- |(k • v)-' [v X (k X Ak)] • 

c J { dp 



(32) 



where the off-resonance interaction with the quasi-static fluctuations has been considered 
(§ 2.2). For the momentum distribution function of (/p) oc |p|~'^' for the statistically 
accelerated particles with E = |p|c (§ 2.4), the RHS of equation (32) is of the order of 
~ [|gA(r)|/(c|p|)] (/p) for j3' ~ 0(1). Therefore, we see that within the present framework, 
the quasi-linear approximation is valid for the test particles with an energy of 3> |q'A(r)|, 
in a confinement region. Note that this relation ensures the condition that the gyroradius for 
the local field strength of |B(r)| greatly exceeds the filament size (coherence length) of order 
~ k~^, namely, E/\qB{r)\ » [equivalently, E ^ |gf74(r)|], except for a marginal region 
of /c ~ R~^. Obviously, this means that in the high-energy regime of £^ ^ |gfA|, the test 
particles are not strongly deflected by a local magnetic fleld accompanying a flne fllament 
with its transverse scale of ~ k^^. On the other hand, in the cold regime of ^ k^li the 
test particles are tightly bound to a local magnetic fleld having the (locally deflned) mean 
strength, violating equation (2). Here it is expected that the bound particles can diffuse 
along the local fleld line, and hence, diffusion theories for a low-/3 plasma are likely to be 
more appropriate, rather than the present theory. 

Summarizing the above discussions, there seem to exist two distinct energy regimes for 
the test particles conflned in the system comprising numerous magnetized fllaments: the 
higher energy regime of E ^ Iq^I, iii which the particles are free from the local magnetic 
traps, and the lower energy regime of E ^ \qA\, in which the particles are bound to the 
local fields, as compared to a low-/3 state. The hierarchy is illustrated in Figure 2, indicating 
the characteristic trajectories of those particles. When shock propagation is allowed, as seen 
in actual AGN jets, the shock accelerator can energize the particles in each energy level. 
At the moment, we are particularly concerned with EHE particle production by a feasible 
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scenario according to which energetic free particles, unbound to small-scale structure of the 
magnetized filaments, are further energized by the shock. In this aspect, the particle escape 
from magnetically bound states, due to another energization mechanism, can be regarded 
as the injection of prcaccelerated particles into the concerned diffusive shock accelerator. If 
the preaccelerator, as well, is of DSA, relying on the gyromotion of bound particles (Drury 
1983; Biermann & Strittmatter 1987; Zank et al. 2004), the preaccelerator also calls for the 
injection (in a conventional sense) in a far lower energy level, owing to, e.g., the Maxwellian 
tail, or the energization of particles up to the energies where the pre-DSA turns on (for 
a review, see Berezinskii et al. 1990). The energy required for this injection, the so-called 
injection energy, could be determined by, e.g., the competition with coUisional resistance. 
In order to distinguish from this commonly used definition of "injection," we refer to the 
corresponding one, owing to the particle escape from the magnetic traps, as the "transition 
injection," in analogy to the bound-free transition in atomic excitation. The energy required 
to accomplish of the transition injection is formally denoted as i^inj ~ |q'^|th, where IgAjth 
represents a threshold potential energy. That is, the particles with q and the energy exceeding 
£^inj are considered to spaciously meander to experience successive small deflection by the 
flelds of many fllaments (Fig. 2), such that the present theory is adequate for describing the 
particle diffusion. This scattering property can be compared to that for the conventional 
QLT in low-/3 regimes: an unperturbed (zeroth order) guiding center trajectory of gyrating 
particles bound to a mean magnetic field must be a good approximation for many coherence 
lengths of particle scatterer. 

If both the injection and transition injection work, the multistep DSA can be realized. 

In the stage oi E <^ \qA{v)\, many acceleration scenarios that have been proposed thus 
far (DSA: e.g., Drury 1983; Biermann & Strittmatter 1987; shock drift acceleration: e.g., 
Webb et al. 1983; some versions of the combined theories: e.g., Jokipii 1987; Ostrowski 
1988; for a review, sec, e.g., Jones & Ellison 1991) can be candidates for the mechanism 
of the preacceleration up to the energy range of £' ~ |gA(r)|, although before achieving 
this energy level, the acceleration, especially for electrons, might be, in some cases, knocked 
down by the energy loss, such as synchrotron cooling, collision with photons, and so on. 
The relevant issues for individual specific situations are somewhat beyond the scope of this 
paper (observability is discussed in § 4). Here we just briefly mention that in the termination 
regions of large-scale jets where the bulk kinetic energy is significantly converted into the 
magnetic and particle energies, a conventional DSA mechanism involving large-scale MHD 
turbulence might work up to EHE ranges (Honda & Honda [2004b] for an updated scenario 
of oblique DSA of protons) . 
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3.2. Timescale of the Diffusive Shock Acceleration 



In the following, we focus on the DSA of energetic free particles after the transition 
injection. Let us consider a typical case of 0i = 0n = 0° in equation (31), reflecting a 
reasonable situation that a shock wave propagates along the jet comprising linear filaments. 
Since the vectors of the random magnetic fields are on the plane transverse to the current 
filaments, this plane is perpendicular to the shock-normal direction. That is, the shock across 
the perpendicular magnetic fields is considered. In this case, no irregularity of magnetic 
surfaces in the shock-normal direction exists, because of k\\ = 0. However, this does not 
mean that the particle flux diffusively across the shock surface is in free-streaming; note 
that the particles crossing the local fields with nonsmall pitch angles suffer the orthogonal 
defiection. 

Anyhow, the injected particles are off- resonantly scattered by the filamentary turbu- 
lence, to diffuse, migrating back and forth many times between the upstream and downstream 
regions of the shock. As a consequence, a small fraction of them can be stochastically ac- 
celerated to very high energy. This scenario is feasible, as long as the filamentary structure 
can exist around the discontinuity, as seen in a kinetic simulation for shock propagation 
(Nishikawa et al. 2003) and an actual filamentary jet (Owen et al. 1989). The timescale of 
this type of DSA is of the order of the cycle time for one back-and-forth divided by the energy 
gain per encounter with the shock (Gaisser 1990). Here the cycle time is related to the mean 
residence time of particles (in regions I and II), which is determined by the diffusive particle 
fiux across the shock, dependent on Kn^^ (eq. [31]). For the moment, the shock speed is 
assumed to be nonrelativistic. Actually, this approximation is reasonable, since the discrete 
knots (for FR I) and hot spots (for FR II), which are associated with shocks (e.g., Biretta et 
al. 1983; Carilli & Barthel 1996), preferentially move at a nonrelativistic speed, slower than 
the speed of the relativistic jets (e.g., Meisenheimer et al. 1989; Biretta et al. 1995). When 
taking the first-order Fermi mechanism into consideration for calculation of the energy gain, 
the mean acceleration time can be expressed as (Lagage & Cesarsky 1983a,b; Drury 1983) 



where Ui and Uu are the flow speed of the upstream and downstream regions in the shock 
rest frame, respectively. The present case of 0^ = (in eq. [31]) provides Kn.c = '^H.C' where 
K\\^i^ is given in equation (26). Here note the relation of Bi — Bu, derived from the condition 
that the current density, ~ V x B^, must be continuous across the shock front. When 
assuming ai — an and /3i = we arrive at the result 




(33) 
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Sv^Trar (r + 1) cE^ 



(34) 



'a,acc — 



8 (a - 1) + 1) {(3a + 2) (r - 1) qlBmU^ ' 



where the definitions oi a = a^, (3a = Pci = -^(i U = Ui — rUu have been introduced. 
Equation (34) is vaUd for arbitrary species of particles "a" having energy Ea, spectral index 
(3a, and charge g^. Note that for the plausible ranges of the values of a, (3a, and r, the value 
of equation (34) does not significantly change. The 0^ dependence is also small, because of 
k ^ 1, refiecting three-dimensional rms defiection of unbound particles (§ 3.1 and Fig. 2). In 
the scaling laws shown below, for convenience we use the typical values of a = 2 (Montgomery 
& Liu 1979) and r — A (for the strong shock limit), although we indicate, in equation (39), 
the parameter dependence of the numerical factor. 



In equation (34) for ions (a ="i"), we set g; = Z\e\ (see footnote 1) and (3i = 3 
(e.g., Stecker & Salamon 1999; de Marco et al. 2003). By balancing equation (34) with 
the timescale of the most severe energy loss process, we derive the maximum possible energy 
defined as £'i,max = Ei. In the environment of astrophysical filaments including extragalactic 
jets, the phenomenological time balance equation can be expressed as 



where tsh, ^i.syn, tn-y, and tnn' stand for the timescales of the shock propagation (§ 3.3.1; 
eq. [36]), the synchrotron loss for ions (§ 3.3.2; eq. [40]), the photodissociation of the nucleus 
(§ 3.3.3; e.g., eq. [43]), and the collision of nucleus "n" with target nucleus "n"' (§ 3.3.4; e.g., 
eq. [46]), respectively. In addition, the energy constraint ascribed to the spatial scale, i.e., 
the quench caused by the particle escape, should also be taken into account (§ 3.3.5). The 
individual cases are investigated below. 



In the actual circumstances of astrophysical jets, the propagation time of a shock through 
the jet, tsh) restricts the maximum possible energy of accelerated particles. The shock prop- 
agation time may be interpreted as the age of knots or hot spots (Honda & Honda 2004b), 
which can be crudely estimated as ~ L/Uprop, where L represents a distance from the central 



3.3. The Highest Energy of an Accelerated Ion 



'i,acc 




(35) 



3.3.1. The Case Limited by the Shock Propagation Time 
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engine to the knot or hot spot being considered and C/prop an average speed of their proper 
motion. When assuming U ~ C/prop, we get the scahng 

U 1 X 10^^-^^ s. (36) 
1 kpc U 

For the case in which the shock is currently ahve as is observed in AGN jets, 4h cannot be 
compared to the "hfetime" of the accelerator that is considered in SNR shocks (e.g., Gaisser 
1990). 

It is mentioned that in AGN jets, the timescale of adiabatic expansion loss might be 
estimated as tad ^ 3L/ (2rUr), where F and Ur represent the Lorentz factor of jet bulk flows 
and the speed of radial expansion, respectively (Miicke et al. 2003). The fact that the jets 
are collimating well with an opening angle of ^ 10° means f/prop/f/r ^ O{10); thereby, 
^sh ^ ^ad ior F < (9(10). Thus, it is sufficient to pay attention to the limit due solely to the 
shock propagation time. These circumstances are also in contrast with those in the SNRs, 
where the flows are radially expanding without coUimation, and the shock propagation time 
(or lifetime) just reflects the timescale of adiabatic expansion loss (e.g., Longair 1992). 

In equation (35), let us ffist consider the case of tj ^cc = ^sh- By equating (34) with (36), 
we obtain the following expression for the maximum possible energy of an accelerated ion: 

^, . 70 f ^) "7^) "7i^) EeV. (37) 

imGVikpcy Vioopcy \0.1cJ ^ ' 

Note the ratio of L/i? ~ 36O/(7r0oa sin 0va) ~ 10 — 100 for the narrow opening angle of AGN 
jets of 0oa ~ 1° — 10° and not-so-small viewing angle 0va (e-g-, for theM87jet, L ~ 23-R— 33i? 
for 0oa ^ 6?9 [Reid et al. 1989] and = 42?5±4?5 [Biretta et al. 1995] or 30° -35° [Bicknell 
& Begelman 1996]). Equation (37) (and eq. [48] shown below) corresponds to the modified 
version of the simple scaling originally proposed by Hillas (1984). 

Concerning the abundance of high-Z elements and their acceleration to EHE regimes, 
the following points 1-4 may be worth noting: 

1. Radial metallicity gradients are expected to be enhanced in elliptical galaxies (e.g., 
Kobayashi 2004). Along with this, a significant increase of heavy elements has been 
discovered in the central region of the nearby giant elliptical galaxy M87 (Gastaldello 
& Molendi 2002), which contains a confirmed jet. 

2. A variety of heavy ions including iron have been detected in a microquasar jet (SS 433; 
Kotani et al. 1996). 
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3. The Haverah Park data favor proton primaries below an energy of ~ 50 EeV, whereas 
they appear to favor a heavier composition above it (Ave et al. 2000) . 

4. The recent Fly's Eye data of ~ 320 EeV are compatible with the assumption of a 
hadron primary between proton and iron nuclei (Risse et al. 2004). 

With reference to this observational evidence, we take the possibility of acceleration of (or 
deceleration by) heavy particles into consideration and indicate the charge (Z) and/or atomic 
number (A) dependence of the maximum possible energies and loss timescales. 



3.3.2. The Case Limited by the Synchrotron Cooling Loss 

The particles deflected by the random magnetic fields tend to emit unpolarized syn- 
chrotron photons, which can be a dominant cooling process. For relativistic ions, the 
timescale can be written as ii^syn — 367r^(A/Z)^[mpC^/(e^S^£'i)], where rup denotes the pro- 
ton rest mass (Gaisscr 1990). In this expression, the energy of an accelerated ion, E^, can 
be evaluated by equating acc with ti^gyn- That is, we have 



ArripC^ 



C{a,pi,r) 



A rUpR fU 



Z"^ ruero \c 



1/3 



(38) 



Here the dimensionless factor ^ is given by 



4^/6(Q;-l) (/5 + !)(/? + 2) (r - 1) 



ar (r + 11 



1/3 



(39) 



and To = e^/ {Anm^c^) stands for the classical radius of the electron, where is the electron 
rest mass. Substituting equation (38) into the expression of ti,syn, the cooling timescale can 
be expressed as a function of the physical parameters of the target object. As a result, we 
find 



syn 



3 X 10 



15 



^_ M y^' ( I mc y ^loopcy/' 



\2Z) 



B 



R 



2/3 



(40) 



Practically, this expression can be used in equation (35) for making a direct comparison with 
the other loss timescales. For example, in the FR sources with B < 1 mG (Owen et al. 1989; 
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Meisenheimer et al. 1989, 1996; Ra^hen & Biermann 1993), we have ti^gyn ^ *sh, so that 
the synchrotron coohng loss is ineffective. It should, however, be noted that in blazars with 
5 > 0.1 G (Kataoka et al. 1999; Miicke & Protheroe 2001; Aharonian 2002), equation (40) 
becomes, in some cases, comparable to equation (36). 

When the equality of ii^acc = ^i,syn is fulfilled in equation (35), equation (38) just provides 
the maximum possible energy of the accelerated ion, which scales as 

^•■--'^"i^n^i^Ti&f ZeV. (41) 

The important point is that ti^^cc and ti^syn are both proportional to B~^, so that the B 
dependence of -Bi^max is canceled out. This property also appears in the case of electron 
acceleration attenuated by the synchrotron cooling (§ 3.4.1). In equation (41), it appears 
that for heavier ions, £^i,max takes a larger value. In the actual situation, however, the 
extremely energetic ions possess a long mfp, and therefore, acceleration may be quenched by 
the particle escape, as discussed in § 3.3.5. 



3.3.3. The Case Limited by the Collision with Photons 



Here we focus on the proton-photon collision that engenders a pion-producing cascade. 
The characteristic time of the collision depends on the target photon spectrum n(eph) in the 
acceleration site, where n(eph) is the number density of photons per unit energy interval for 
photon energy Cph- For n(eph) oc (e.g., Bezler et al. 1984), typical for the FR sources 
(Rachen & Biermann 1993), the timescale can be expressed as tp^ ~ [wm/(xWph)]tp,syni where 
X ~ 200 for the average cross section of cr^p ~ 900 ^barns (Biermann & Strittmatter 1987), 
Mph denotes the average energy density of target photons, and tp,syn = ^i,synU=z=i (the 
subscript "p" indicates proton). Thus, the expression of tp^ includes Ep, i.e., the energy of 



the accelerated proton. This can be evaluated by equating tp^a 
of 



with tpy, to have the form 



Ep e(«,A,r) 



nip R fU 
me ro 



1/3 



(42) 



where the definition ry,, = Uph/Mm has been introduced. Substituting equation (42) into the 
expression of tp^, we obtain the following scaling of the photomeson cooling time: 
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Note that for rju — x"^ ~ 5 x 10"^, we have tp^ = tp,syn- 

If the equahty of tp^acc = is satisfied in equation (35), then equation (42) gives the 
maximum possible energy, which scales as 

For rju — x"^: equation (44) is identical with equation (41) for ^4 = Z = 1. 



3.3.4- The Case Limited by the Collision with Particles 

The nucleus-nucleus collisions involving spallation reactions can also be a competitive 
process in high-density regions. For proton-proton collision, the timescale can be simply 
evaluated by tpp' — (rip/CTpp/c)"^, where Up' is the number density of target protons, and 
CTpp/ ~ 40 mbarns denotes the cross section in high-energy regimes. The timescale can be 
rewritten as 

1 CII1~^ 

tpp> ~ 8.3 X 10^^ s. (45) 

rip' 

It is found that for tenuous jets with rip/ <S 1 cm~^, the value of equation (45) is larger than 
the conceivable value of equation (36); that is, the coUisional loss is ineffective. 

For the collision of an accelerated proton with a nonproton nucleus, the timescale can be 
evaluated by the analogous notation, ipN' = (n^'CpA'c)"^, where ua' is the fractional number 
density of the target nuclei having atomic number A' > 1. Here we use an empirical scaling 
of the cross section, apA' ~ nr^A''^^^, where Tq ~ 1.4 x 10~^^ cm, although the value of Tq 
may be an overestimate for very high energy collisions (e.g.. Burbidge 1956). Combining 
tpN' with tpp/, in general the timescale for collision of a proton with a nucleus of an arbitrary 
composition can be expressed as 
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where Up' and ua' are both in units of cm ^. 

In equation (35), we consider the case of tp,acc — tpn'- By equating (34) with (46), we 
obtain the following expression for the maximum possible energy of an accelerated proton: 



/ 100 cm-:* B f R t/ 

As for the collision of an arbitrary accelerated nucleus with a target nucleus, we can anal- 
ogously estimate tnn' and -Ei^max- In particular, the heavier nucleus-proton collision is more 
important, since its timescale tnp' is of the order of tpp//A^/^: for larger A and rip/, it can 
be comparable to the other loss timescales. For example, the parameters of A = 56 (iron) 
and rip' ~ 100 cm~^ lead to inp' ~ 4 x 10^^ s. For the case of fi^acc = ^np' in equation (35), 
we have the scaling of E'i^max ~ 0.6Z^/^(2Z/74)^/^£'p^inax, where -Bp^max is of equation (47) for 
V >EA'>l^A'>l"/^■ 



5.5.5. Quenching by Particle Escape 

The particle escape also hmits its acceleration; that is, the spatioscale of the system 
brings on another energy constraint. Relating to this point, in § 2.4 we found the relation of 
K ~ 1, meaning that the anisotropy of the spatial diffusion coefficient is small. It follows that 
the radial size of the jet (rather than L) affects the particle confinement. Recall here that 
in the interior of a jet the magnetic field vectors tend to be canceled out, whereas around 
the envelope the uncanceled, large-scale ordered field can appear (Honda & Honda 2004a). 
From the projected view of the jet, on both sides of the envelope the magnetic polarities are 
reversed. 

The spatially decaying properties of such an envelope field in the external tenuous 
medium or vacuum might influence the transverse diffusion of particles. The key property 
that should be recalled is that for r ^ distant from a filament, the magnetic field 
strength is likely to slowly decay, being proportional to ~ {kr)^^ (Honda 2000; Honda & 
Honda 2002). It is, therefore, expected that as long as the radial size of the largest filament, 
i.e., correlation length, is comparable to the radius of the jet (§ 2.3), the long-range field 
pervades the exterior of the jet, estabhshing the "magnetotail" with the decay property of 
~ (^min'^)~^ for r » A;~?^(~ R). In fact, in a nearby radio galaxy, the central kiloparsec-scale 
"hole" of the inner radio lobe containing a jet is filled with an ordered, not-so-weak (rather 
strong) magnetic field of the order of 10 — 100 /iG (Owen et al. 1990), whose magnitude 
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is comparable to (or ~ 10 % of) that in the jet (Owen et al. 1989; Heinz & Begelman 
1997). Presumably, the exuding magnetic field plays an additional role in confining the 
leaky energetic particles with their long mfp of Ax(~ c^^l'^^-i) ~ -R- 

In this aspect, let us express an effective confinement radius as — pR, where p ^ 1, 
and impose the condition that the accelerator operates for the particles with the transverse 
mfp of \±< Rc- Then the equality gives the maximum possible energy in the form of 



E,,^^^2m Z~py \ ^ J EeV. (48) 
1 mG 100 pc 

Values of E'i^max (and £'p,max) derived from the time balance equation (35) cannot exceed 
that of equation (48). It appears that equation (48) can be compared to the energy scaling 
derived from, in the simplest model, the maximum gyroradius in a uniform magnetic field. 



3.4. The Highest Energy of an Accelerated Electron 

In a manner simiar to that explained in § 3.3, we find the generic scaling for the achiev- 
able highest energy of electrons. In equation (34) for electrons (a ="e"), we set — — |e| 
and Pe = '2 (e.g., Meisenheimer et al. 1989; Rachen & Biermann 1993; Wilson & Yang 2002). 
By balancing equation (34) with the timescale of a competitive energy loss process, we derive 
the maximum possible energy defined as £^e,max = Eg. The time balance equation can be 
written as 



^e,acc — min (ie,syn) ^io ^br) ) (^9) 

where te.syn, ^ic, and tbr stand for the timescales of the synchrotron loss for electrons (§ 3.4.1; 
eq. [51]), the inverse Compton scattering (§ 3.4.2; eq. [54]), and the bremsstrahlung emission 
loss (§ 3.4.3; eq. [57]), respectively. For positrons the method is so analogous that we omit 
the explanation. 



3.4- 1- The Case Limited by the Synchrotron Cooling Loss 

For electrons, the synchrotron cooling is a familiar loss process, and the timescale can 
be expressed as te,syn — SGir'^mjc'^ /(c'^B'^Eq). In this expression, the energy of an accelerated 
electron, E^, can be evaluated by equating ie,acc with ^e.syn, to give 
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R ru 

ro V c 



1/3 



(50) 



Substituting equation (50) into the aforementioned expression of tc,syn) the coohng timescale 
can be expressed as a function of the physical parameters of the target object: 



B 



R 



U 



(51) 



This can be used in equation (49) for comparison with the other loss timescales. 



When the equality of ie.acc = ^e,syn is satisfied in equation (49), equation (50) gives the 
maximum possible energy, which scales as 



Ee,max ~ 50 



R 



100 pc 



1/3 



U 



2/3 



PeV. 



(52) 



According to the explanation given in § 3.3.2, equation (52) is independent of B (see also 
eq. [41]). The striking thing is that for plausible parameters, the value of -Ee,max is significantly 
larger than that obtained in the context of the simplistic QLT invoking the Alfven waves 
(Biermann & Strittmatter 1987). This enhancement is, as seen in equation (28), attributed to 
the smaller value of the diffusion coefficient for electrons, which leads to a shorter acceleration 
time, i.e., a smaller value of equation (33), and thereby to a higher acceleration efficiency. 



3.4-2. The Case Limited by the Inverse Compton Scattering 

For the case of > 1, the inverse Compton scattering of accelerated electrons with tar- 
get photons can be a dominant loss process. Actually, the environments of AGN jets often 
allow the synchrotron self-Compton (SSC) and/or external Compton (EC) processes. The 
characteristic time of the inverse Comptonization can be estimated as tic ~ (^e,syn/?7u)(cT/cKN), 
where ctt = 87rrQ/3 and (TKN(eph-E'e) denote the total cross sections in the Thomson limit of 
CphE'c ^ m^c^ and the Klein-Nishina regime of Cph-Ec ^ m^c^, respectively (e.g., Longair 
1992). The expression of tic includes E^, which is determined by numerically solving the 
balance equation of te,acc = ^ic- Because of ctkn < ct, the value of E^^ is found to be in the 
region of 
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R fU^ ^ 



(53) 



in the whole range of Cph- Note that the equaUty in equation (53) reflects the Thomson limit 
of ctkn/ct = 1- Substituting the value of E^ into the expression of tic, we can evaluate the 
scattering time, which takes a value in the range of 

> 5 X io»,y3lO-'°argcm-3 / lOOpcX /0_lcX ^ ^^^^ 

■Uph \ R J \ U J 

For a given parameter the larger value of ?7,ph depresses the lower bound of tic, though 
the Klein-Nishina effects prolong the timcscalc. It should be noted that the evaluation of tic 
along equation (54) is, in equation (49), meaningful only for > 1; that is, the relation of 
riu<l ensures tic > te,syn- 

For the case of tc,acc = tic in equation (49), -Ee, conforming to equation (53), gives the 
maximum possible energy, which takes the value of 

/ 1 \ 1/3 / p \^/^ f TJ \ 

£;e,.a.>50 - -— — - PeV (55) 



for ?7m > 1. Again, note that the Thomson limit sets the lower bound of £'e,max- It is found 
that the Klein-Nishina effects enhance the value of -Ecmax in the regime of Cph ^ m^c"'/-E'c,max- 
Note here that -Ec,max cannot unlimitedly increase in actual circumstances but tends to be 
limited by the synchrotron cooling. Combining equation (52) with equation (55), therefore, 
we can express the allowed domain of the variables as follows: 

Vu 



Note that the upper bound reflects the synchrotron limit. In the critical case of = 1 
reflecting the energy equipartition, the generic equation (56) degenerates into equation (52). 



3.4-3. The Bremsstrahlung Loss 



The bremsstrahlung emission of electrons in the Coulomb fleld of nuclei whose charge 
is incompletely screened also affects the acceleration. The timescale can be evaluated by 
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the notation t^r — (jT-z'Crad.eZ'c)"^, where nz' is the fractional number density of the target 
nuclei having charge number Z' and (Jrad.eZ' describes the radiation cross section (e.g., Heitler 
1954). When the screening effects are small, for interaction with a heavy composite we have 



where nz' is in units of cm ^. 

In the peculiar environments of high density, enhanced metallicity, and lower magnetic 
and photon energy densities, equation (57) may be comparable with equation (51) or (54). 
In ordinary AGN jets, however, the corresponding physical parameters seem to be marginal. 
Note also that the bremsstrahlung timescale for ion-ion interactions is larger, by the order of 
{A^ /Z'^){mp/'mey ~ 10^/Z^, than the value of equation (57), and found to largely exceed the 
value of equation (36), namely, the age of the accelerator. That is why the ion bremsstrahlung 
has been excluded in equation (35). 



The feasibility of the present model could be verified by the measurement of energetic 
photons emanating from a source, typically, bright knots in nearby AGN jets. In any case, the 
electrons with energy -Ee.max, given in equation (56), emit the most energetic synchrotron 
photons, whose frequency may be estimated as u* ~ {Ef.^^a_^/meC^Y{eB /rriQc), where the 
mean field strength B has been compared to the rms strength B. For Ee,max ~ 10 PeV as 
an example, the frequencies of u* > 10^^ Hz are found to be achieved for B > 10 //G. In 
the gamma-ray bands, however, the energy fiux of photons from the synchrotron originator 
is predicted to be often overcome by that produced by the inverse Comptonization of target 
photons. In this case, as far as the condition of eph-Ec,max ^ {fncC?Y is satisfied, the boosted 
photon energy is given by epj^ ~ -Ec,max, independent of the target photon energy Cph, thereby 
irrespective of SSC or ECs. This is in contrast to another case of eph-E'e,max ^ (^eC^)^) in 
which one has e'^-^ ~ eph(ii^e,max/*^eC^)^ dependent on eph. Apparently, for the extremely high 
energy ranges of i^e.max achieved in the present scheme, the former condition is more likely 
satisfied for a wide range of eph. Therefore, in the circumstances that the source is nearby such 
that collision with the cosmic infrared background, involving photon-photon pair creation, 
is insignificant, e'^-^ (~ -E'e,max) just gives the theoretical maximum of gamma-ray energy, 
although the Klein-Nishina effects also take part in lowering the flux level. This means, in 
turn, that a comparison of the e'j^ value (multiplied by an appropriate Doppler factor) with 




(57) 



4. 
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the observed highest energy of the Compton emissions might constitute a method to verify 
the present DSA for electrons. 

The case for this method is certainly solidified when the operation of the transition 
injection (§ 3.1) is confirmed. Making use of the inherent property that the synchrotron 
photons emitted by electrons having an energy above |e^(r)| reduce their polarization, the 
energy hierarchy can be revealed by the polarization measurements, particularly, with wide 
frequency ranges and high spatioresolution. According to the reasoning that the critical 
frequency above which the measured polarization decreases, t'c(i"), ought to be of the or- 
der of ~ [|e/l(r)|/mcC^]^[|eB(r)|/mcc], the related coherence length can be estimated as 

~ c{z/c(r)[mec/|eB(r)|]'^}-'^/^. Note that when the locally defined gyroradius reaches this 
critical scale, the bound electrons are released. In actual circumstances, u^. and the polariza- 
tion for a fixed frequency band are, if anything, likely to increase near the jet surface, where 
the large-scale coherency could appear (§ 3.3.5). This may be responsible for the results of 
the polarization measurement of a nearby filamentary jet, which indicate a similar transverse 
dependence (Capetti et al. 1997). In the sense of £'e,max/£'mjl9=-|e| ^ -Si,max/£'inj|g=z|e|) the 
transition injection condition for electrons is more restrictive than that for ions. Thus, ob- 
servational evidence of the present DSA scenario for energetic electrons will, if it is obtained, 
strongly suggest that the same scenario operates for ion acceleration, providing its finite 
abundance. 

To summarize, we have accomplished the modeling of the diffusive shock accelerator 
accompanied by the quasi-static, magnetized filamentary turbulence that could be self- 
organized via the current filamentation instability. The new theory of particle diffusion 
relies on the following conditions analogous to those for the conventional QLT: (1) the test 
particles must not be strongly deflected by a flne fllament but suffer the cumulative small 
deflection by many fllaments, and (2) the transverse fllament size, i.e., the coherence length 
of the scatterer, is limited by the system size transverse to the fllaments; whereas, more 
importantly, it is dependent on neither the gyration, the resonant scattering, nor the explicit 
limit of the weak turbulence. We have derived the diffusion coefficient from the quasi-linear 
type equation and installed it in a DSA model that involves particle injection associated 
with the bound-free transition in the fiuctuating vector potential. By systematically taking 
the conceivable energy restrictions into account, some generic scalings of the maximum en- 
ergy of particles have been presented. The results indicate that the shock in kiloparsec-scale 
jets could accelerate a proton and heavy nucleus to 10 — 100 EeV and ZeV ranges, respec- 
tively. In particular, for high-Z particles, and electrons as well, the acceleration efficiency 
is significantly higher than that derived from a simplistic QLT-based DSA, as is deduced 
from equation (28). Consequently, the powerful electron acceleration to PeV ranges becomes 
possible for the plausible parameters. 
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We expect that the present theory can be, mutatis mutandis, apphed for solving the 
problem of particle transport and acceleration in GRBs (Nishikawa et al. 2003; Silva et al. 
2003). The topic is of a cross-disciplinary field closely relevant to astrophysics, high-energy 
physics, and plasma physics involving fusion science; particularly, the magnetoelectrodynam- 
ics of filamentary turbulence is subject to the complexity of "flowing plasma." In perspective, 
further theoretical details might be resolved, in part, by a fully kinetic approach allowing 
multiple dimensions, which goes far beyond the MHD context. 



A. CALCULATION OF THE INTEGRAL COMPONENTS 7i2, hi, AND h2 

For instruction, we write down the derivation of equations (24) and (25) for the col- 
lisionless scattering of injected test particles by magnetized current filaments having the 
configuration illustrated in Figure 1. Making use of equation (10), hj for ij ^ 11 (eq. [12]) 
can be explicitly written as 

= / d^kdcl (k . V) A I _ (k . v)]-^ (k . ^) (/p) } , (Al) 

hi = j d'\^MK,.v\\ (k • ^) |[(^ - (k • v)]-^ (k • v) ^ (/p)} , (A2) 
/22 = ^^ / d^kd^|^|^,,.|| (k-^) {[^-(k-v)]-%|| (k-^) (/p)}, (A3) 

where we have used a standard correlation function (eq. [13]) reflecting random magnetic 
fluctuations on the transverse plane to the linear current fllaments (see Fig. 1). Recalling 
the causality principle and noticing that the real part does not contribute to the integration, 
we get 

/i2 = / ^'^MA\l^ • 4 ~ ■ ■ ^) ' ^^^^ 

hi = / d^kd^l^l^,,.,, (k • ^) 1^ [^ - (k • v)] (k • v) ^ (/p)} , (A5) 
/22 = ^ y d^kd^|^|^,,.|| (k • 1^ [^ - (k • v)] (k • ^) (/p)} . (A6) 
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Again, we use an ad hoc equation (16), valid for a quasi-static mode that retains the narrow 
spectral peak around a; ~ 0. Assuming that the magnetic turbulence is isotropic on the 
transverse plane, the angular average of equations (A4)-(A6) is carried out. Taking account 
of the off-resonant scattering of particles by the quasi-static random fields gives 



where the ordering and the definition of Ik are the same as those denoted in § 2.2. For 
a given momentum distribution function of (/p) oc IpI"^'^^^^ for the test particles with an 
ultrarelativistic energy of = |p|c = c(py the partial derivatives can be estimated 

as 



dp\\ \vx_ 



dp 



2 (/p) , (AlO) 



(/p) ~ (/3 + 2)(/3 + 4)V'iV'2|^(/p> 



dp±dp\\ 

d ( v\\ \ d 



dp± \v± J dp 



(/p) , (All) 



where V'l = P||/|p| and ■02 = P±/|p|- Using equations (AlO) and (All), and equation (22) 
concerning the /c-spacc integration for k^gx 00, one can arrange equations (A7)-(A9) in 
the form of 1^ ~ (/p) to obtain the expressions of equations (24) and (25). 
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Fig. 1. — Schematics of the midscale configuration of zeroth-order current density vectors 
J ~ Jz {dotted arrows along the "poles" representing current filaments) and magnetic field 
vectors B = (5^,., 5^, 0) {bold solid arrows randomly distributed on the shaded plane), and 
a spatial profile of the vector potential A = Az {bold curve), embodying the magnetized 
current filaments that constitute the bulk of a large-scale jet. Here the scalar A{x,y) has 
been depicted as a function of x for a fixed y. [See the electronic edition of the Journal for 
a color version of this figure.] 
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Fig. 2. — Schematics of the trajectories of a test free particle with energy exceeding a 
threshold potential energy |gA|th, and a gyrating, bound particle with energy below the local 
potential \qA{x,y)\ [light solid curves with arrow). In the energy hierarchy, the transition 
from the bound state to the free state can be compared to the particle injection for the 
present DSA, which requires the unbound particles to suffer successive small deflections by 
the random magnetic fields. For an explanation, see the text, [^ee the electronic edition of 
the Journal for a color version of this figure] 



